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ABSTRACT 

I We present the first automated spectroscopic search for disk-galaxy lenses, using the 

^ I Sloan Digital Sky Survey database. We follow up eight gravitational lens candidates, 

■ selected among a sample of ~ 40000 candidate massive disk galaxies, using a combina- 

tion of ground-based imaging and long-slit spectroscopy. We confirm two gravitational 
lens systems: one probable disk galaxy, and one probable SO galaxy. The remaining 
systems are four promising disk-galaxy lens candidates, as well as two probable grav- 
itational lenses whose lens galaxy might be an SO galaxy. The redshifts of the lenses 
are ziens ~ 0.1. The redshift range of the background sources is ^source ~ 0.3 — 0.7. 
The systems presented here are (confirmed or candidate) galaxy-galaxy lensing sys- 
^ ' tems, that is, systems where the multiple images are faint and extended, allowing an 

^ ■ accurate determination of the lens galaxy mass and light distributions without contam- 

ination from the background galaxy. Moreover, the low redshift of the (confirmed or 
. candidates) lens galaxies is favorable for measuring rotation points to complement the 

i lensing study. We estimate the rest-frame total mass-to-light ratio within the Einstein 

QQ I radius for the two confirmed lenses: we find Mtot/Lj = 5.4 it 1.5 within 3.9 it 0.9 kpc 

^ ; for SDSS J081230.30+543650.9, and Mtot/Lj = 1.5 ±0.9 within 1.4 ±0.8 kpc for SDSS 

J145543.55±530441.2 (all in solar units). Hubble Space Telescope or Adaptive Optics 



o 



o 

oo 
O 



imaging is needed to further study the systems. 

Subject headings: dark matter — galaxy: fundamental parameters (M/L) — galaxy: 
spiral — galaxy: structure — gravitational lensing 



1. Introduction 

The Lambda-Cold Dark Matter (ACDM) model has gained the place of cosmological paradigm 
to describe our Universe, yet its success at explaining the large-scale universe is not reproduced 
at galactic scales. The central cusps of simulate d dark matter (DM) halos are in disa greement 



with the observed cores of spiral galaxy DM halos (jGentile et al.ll2004l : iKassin et al.ll2006l ). and the 



disk galaxies formed in ACDM cosmological simulations appear to be offset from the Tully-Fisher 
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relation (TF; Tullv fc Fisheij 1977 ) observed for distant disk galaxies ( Navarro Steinmetz 2000al lbl: 
Portinari &: Sommer-LarsenI 120071 : iDutton et al.ll2007l ). Gravitational lensing provides a promising 
way of constraining the mass distribution of disk galaxies and measuring their mass-to-light (M/L) 
raticQ. These are key components for comparing the mass of disks grown in cosmological simulations 
to the luminosity of observed galaxies. 

It is still unclea r wheth er the TF relation zer o point problem originates in the simulations 
(IDutton et alj (|2007l ): but see lGovernato et alj (|2007l )). in the ACDM cosmology (jvan den Bosch et al, 
20031 ). or in the generaliza.t ion of the Milky Way (MW) initial mass function (IMF) to other disk 
galaxies (jFlvnn et al.ll2006l l. However, recent studies found the MW to be offset from the observed 
TF relation by about the same amount as disk-galaxi es simulated using Solar Neighborhood IMFs 
(jFlynn et al.ll2006l : iPortinari &: Sommer-LarsenI 120071 ). This suggests that the simulations may be 



consistent with t he model they r e 



disk galaxies (see IPortinari et al.l (|2004l ) for a review of IMFs) . 



produce, but the Galactic IMF may not be representative of all 



The M/L ratio of disk galaxies can be used to cons train the shape of the IMF (de Jong &: Bell 



20071 ). as well as to test the maximum disk hypothesis (jvan Albada &: Sancisilll986l ). Yet, to date, 
detail ed individua. 1 meas urements of M / L ratios exist for three disk galaxies only. These being the 



MW ( Flvnn 



et al 



(I2OO6I) . Mj./L f = 1.20 for the stellar matter in the Galactic disk), the Sc galaxy 
NGC 4 414 (IValleio et al.l hmi \ M^/Lj ~ 1), and the Sab s piral lens 22 3 7-H03 05. i.e., the Einstein 
Cross dTrott fc Websterl (jiooi), M^/Lj = 1.1 for the disk: iTrott et al.l ^QM), M^/Lb = 1.2 for 



the disk). There are several methods used to determine indirectly the M/L ratio of disk galaxies, 
for example, stellar population models (assuming the shape of the IMF), and relati ons between the 
color and M/L ratio der ived from the maximum rotation curves of spiral galaxies (jBell &: de Jong 



2001 



Salucci et al.ll2008l ^. 



Gravitational lensing can provide an independent measurement of the M/L ratio of disk galax- 
ies. The lensing geometry is sensitive to the total mass enclosed within the Einstein radius, and 
to the total projected ellipticity, giving information on the combination of the bulge, disk and DM 
halo. However, a disk-halo degeneracy r emains due to the lack o f co nstraints on the hal o ellipticity. 
The study of the spiral lens B1600-F434 (jjaimsen fc Hiorthlll997l ;i bv lMaller etTaP (I2OO0I ) proved the 
feasibility of breaking this degeneracy assuming an independent mass measurei nent which would 
constr ain t he disk and h alo co ntributions at larger radii. F ollowing this metho d, iTrott &: Webster 
(j2002l ) and iTrott et al.l (j2008l ) studied the Einstein Cross (jHuchra et al.l Il985l ) and measured the 
bulge and disk M/L ratios of the spiral lens, using in addition to the lensing constraints two H I 
rotation points, as well as rotation curve and velocity dispersion profiles. These th ree studies found 
evidence for a sub- maximum disli^, in contradiction with rotation curve studies (Salucci &: Persic 



19991 ). 



^The M/L ratios are expressed in solar units throughout the paper. 

^The theoretical study of Shin fc Evans! 1 2007 1 also points towards spiral lens galaxies having sub-maximum disks. 
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Unfortunately, only ten confirmed disk-galaxy lenses, and three late-type galaxy lenses which 
might have a disk, are known to date (see Table [1]). In comparison, about a hundred elliptical- 
galaxy lenses we re already knowrl^ before the large surveys of recent years ( Bolton et al. 2008 : 
Faure et al.ll2008l ). This difference is due in part to the fact that disk galaxies are less massive than 
elliptical galaxies, and so have a lower multiple image cr oss-section. Theoretical models predict 



that only 10% to 20% of lenses are due to spiral galaxies (jKeeton &: Kochaneklll998l : iMoller et al 



20071 ). Moreover, their detection by optical imaging can be made more challenging due to the small 
separation of the lensed images and by dust extinction in the disk. 

Fortunately, spectroscopic selection of candidates opens new perspectives for finding disk- 
galaxy lens systems. The Sloan Lens Advanced Camera for S urveys (SLACS) pro j ect h as pioneered 



this new technique to discover strong gravitational lenses (iBolton et al 



2004, 



l200fil^. u sing the 

large spectroscopic database of the Sloan Digital Sky Survey (SDSS; lYork et al.l l200d ). They 



selected luminous red galaxies with absorption-line dominated spectra that also showed at least 
three emission lines from a background galaxy along the line-of-sight. Imaging with the Hubble 
Space Telescope (HST) has confirmed over 70 new gravitational lensejfl, with a survey efficiency 
over 65 %. All of the confirmed systei ns had measured le ns and source redshifts, which is required 
for dete rmining the mass of the lens ([Bolton et al.l 120081 ). The Optimal Line-of-Sight (OLS) lens 



survey (jWillis et al. 



2005 



20061 ) extended the search to those lens candidates showing only one 



emission line from a background galaxy. This allowed lensed galaxies at higher redshifts to be 
investigated. They found seven new early- type gravitational lenses from the SDSS, which they 
confirmed using a combination of ground-based imaging and spectroscopic observations. 

The field of strong-lens searching has experienced a dramati c surge in recent yea,rs, from the 
discovery of new radio- loud or optical- loud lensed quasars (^e.g.. iBrowne et al.l 20031: iMvers et al 



2003 



2008a 



Oguri et al. 



2006 



Kavo et al 



2007 



Ghosh &: Narasimha 



2008; 



Inada et al 



2008: 



b|), to seredenpitous galaxy-galaxy lens discoveries f e . g. . Fassnacht et al. 2006; B elokurov et al 



2007 ) , strong lensing opti c al and spectroscopic surveys (e.g.. IWillis et al.ll2005i . .2006. : .Bolton Sz Buries 



Oguri et al 



20071 : iBolton et al.l l2008l : iFaure et al.l l2008l : iLimousin et al.l l2008l^. as well as the study ofthe 
methodology and biases of strong-lens automatic searches (e.g., iMoller et al.l 20071: lOguril 120071 : 
Dobler et al]l2008l : iMandelbaum et al.l[20(M : Marshall et al.lboosi : iNewton et al]l2008l ). 



Spectroscopic selection can be used to efficiently find disk-galaxy lenses in the SDSS, building 
on techniques developed by the SLAGS and OLS-Lens surveys. However, the foreground galaxy as 
well as the background galaxy can have strong emission lines, making the search more difficult. Also, 
the small number of disk-galaxies discovered to date indicates that this project can be challenging. 
Yet, here we show that it is feasible. 



^CASTLES, |http : //cf a-www . harvard . edu/castlesT] 

^We note that the full SLACS sample of confirmed lenses presented in Bolton et al. 1 20081 ) contains six systems 
classified as late-type (Sa or later). 
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In addition, galaxy-galaxy lenses bring better constraints than lensed quasars on the mass 
distribution of the lens galaxy. The lensed images are typically extended, containing much more 
information than point sources, and they are faint, allowing accurate photometry and astrometry 
of the disk galaxy to be carried out. Moreover, the low redshift of z ~ 0.1 for SDSS galaxies is 
favorable for adding rotation curve measurements to the lensing constraints in the future. 

The purpose of this paper is to prove the feasibility of efficiently discovering disk-galaxy lenses 
by an automated spectroscopic search, using the large spectroscopic database of the SDSS. This is 
not a statistical study. We present a method to find low redshift galaxy-galaxy lenses where the 
lens galaxy is a disk, these systems being particularly well suited for measuring the M/L ratios of 
disk galaxies and the mass distribution in their central parts. 

The paper is organized as follows. The selection of the lens candidates is explained in § [2j 
In § [3] we present the ground-based follow-up of the lens candidates with imaging and long-slit 
spectroscopy; with the combined results we identify the genuine disk-galaxy lenses. In § d] we 
estimate the M/L ratio of the confirmed lenses. Improvements and perspectives are discussed in 
§ O Conclusions follow in § [6l 

Throughout the paper we assume a ACDM cosmology with the following cosmological param- 
eter values: Q.m = 0.3, = 0.7, Hq = 70 km s~^ Mpc~^. In this cosmology, 1" corresponds to 
1.84 kpc at a redshift z = 0.1. 



2. Lens candidates selection 



2.1. Massive disk galaxies in the SDSS 



The SDSS provides a large optical ima ging and spectros copic database of astronomical objects 



covering more than a quarter of the sky ( York et al. 200o l) . Defining proper selection criteria to 



constitute a sample of massive disk galaxies out of the SDSS database is not trivial. Disk galaxies 
have very diverse spectral types, ranging from the absorption dominated spectra of SO and Sa 
galaxies to the emission line spectra of S b and Sc galaxies. Ho wever, it is possible to identify disk 
galaxies using color criteria. The study of lStrateva et al.l (120011 ) shows that galaxies in the SDSS are 
distributed in two groups divided by the color separator u — r = 2.22 around redshift z = 0.1, with 
a redshift evolution roughly parallel to the separator line. For each population, the peak density 
corresponds to elliptical and spiral galaxies, respectively. This provides a good criterion to select a 
sample of spiral galaxies in the SDSS. However, because gravitational lensing depends highly on the 
mass of the lensing galaxy, and spiral galaxies are on average less massive than elliptical galaxies, 
we want to find massive disk galaxies to increase the probability that the background galaxies will 
be multiply imaged. That is, try to select preferentially SO and Sa galaxies, which are more massive 
than Sb and Sc galaxies. 



The study of 



Fukugita et al 



(|l995l ) gives more constraints on the color of different galaxy types 
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in the SDSS passband system. Tables are available for various galaxy redsliifts {z = 0.0, 0.2, 0.5, 0.8); 
we use the z = 0.2 galaxy color values, as we impose on our sample the redshift range 0.1 < z < 0.3. 
The upper redshift l imit is fixed by the redshift selection in the Main Galaxy Sample of the SDSS 
([Percival et al.l 120071 ). The lower redshift limit is chosen to increase the efficiency of our survey. 
Indeed, about half of the SDSS Main Galaxy Sample lies below z = 0.1, but these galaxies have 
a lower pro bability of being grav itational lenses due to their low redshift. We cor npare the color 



values from lFukugita et al.l ( 19951 ) to the color distribution of SDSS galaxies found by lStrateva et al 



(2001 



, see their Fi g . 1). The color of galaxies at z 



= 0.2 calculated in the SDSS passband system 
by iFukugita et al.l (|l995l ) have a higher (redder) u — q when compared to those galaxies actually 
observed by the SDSS in the same passband system (IStrateva et al.ll200ll). Th erefore, we decided 



to use the color of Sab and Sbc galaxies at z = 0.2 from IFukugita et al.l (jl995l ). which cover part 
of the region of color-spac e situated between t he density peaks for the elliptical and spiral galaxy 
populations in the SDSS (jStrateva et al.ll200ll ). This is where we expect to find the most massive 



disk galaxies in color-space. 

This leads us to define a color selection criteriorH 1.4 < u—g < 2 and 0.7 < g—r < 1, in addition 
to the r < 17.7 selection intrinsic to the SDSS Main Galaxy Sampl e. We find a sample of 41625 can - 
didate massive disk galaxies, using the SDSS Fifth Data Release (lAdelman-McCarthv et al.l 120071 ). 
Being based only on color, this selection will include some elliptical galaxies. Visual inspection 
and galaxy modeling will be necessary to identify the disk galaxies among the gravitational lens 
candidates. 



2.2. Spectroscopic selection 



The selection of lens candidates among our sample of massive disk galaxies is two-stepped. 
First, we select spectra presenting evidence for at least three emission lines from a background 
galaxy. Second, we filter these candidates based on redshift and magnitude criteria to keep only 
those candidates with the highest lensing probabilities. 

We developed an automated method to select the spectra presenting evidence of at least three 
emission lin es from a backgroun d galaxy. The method used here is similar to that u sed in the SLAGS 
project (see iBolton et al.l l2004l ) and the OLS-Lens survey (see lWillis et al.ll2005l ). but it differs in 
the details. While these two works performed independent processing and b ackground subtrac tion 
on the SDSS spectra (a detailed comparison of the methods is presented in lWillis et al.ll2005l ). we 
used directly the processed spectra distributed by the SDSS, particularly the continuum-subtracted 
and noise spectra. The main features of the continuum-subtracted spectra are skyline residuals left 
from skyline subtraction, galaxy emission lines (if the foreground galaxy has an emission spectrum). 



^We used a combination oi u — g and g ~ r color rather than u — r color in order to match the search criteria 
available in the SDSS SkyServer Spectroscopic Query Form, which we used to retrieve and download our sample of 
galaxies. 



Table 1. Disk-galaxy lenses 



Name 




-^Icns 


^source 


'^Einstein 


disk M/L 


Comments 


ref 


B0218+357 


2 


0.685 


0.944 


0.17 




Too crowded by the quasar images 


1, 2 


B1600+434 


2 


0.41 


1.59 


0.70 






3, 4, 5 


CXOCY J220132. 8-320144 


2 


0.32 


3.90 


0.41 


KU/Lv = 4 


One [O II] rotation point 


6 


OAC-GL J1223-1239 


2 


0.4656 


unknown 


0.42 






7 


PKS 1830-211 


2 


0.886 


2.507 


0.491 




System near to galactic plane 


8 


PMN J2004-1349 


2 


unknown 


unknown 


0.56 






9 


Q2237-(-0305 


4 


0.0394 


1.695 


0.9 


A/./Lj = 1.1 
mJ/Lb = 1.2 


Two H I rotation points and 
stellar-kinematic information 


10, 11, 12 


SDSS J0841-I-3824 




0.1159 


0.6567 


4.21 






13 


SDSS J1432-I-6317 




0.1230 


0.6643 


5.85 






13 


SDSS J2141-0001 




0.1380 


0.7127 


1.81 






13^= 



''Number of lensed images. 
''In arcsec. 

'^Three other late-type galaxy lenses are published in Bolton et al. (2008), but whether the lens galaxy is a disk or an irregular galaxy is not 
confirmed as yet. 

References. — (1) Patnaik et al. 1992; (2) York et al. 2005; (3) Jaunsen & Hjorth 1997; (4) Koopmans et al. 1998; (5) Mailer et al. 2000; 
(6) Castander et al. 2006; (7) Covone et al. 2008; (8) Winn et al. 2002; (9) Winn et al. 2003; (10) Huchra et al. 1985; (11) Trott & Webster 
2002; (12) Trott et al. 2008; (13) Bolton et al. 2008 
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Fig. 1. — Distribution of the foreground galaxy redshift zp, the background galaxy redshift zb 
and the r magnitude of our eight final lens can didates compared to the SLACS confirmed strong 
lenses. The sample of SLACS strong lenses from lBolton et al.l (j2006l ) was the one we actually used 
to select our fin a l lens candidates. We show also the final sample of SLACS strong lenses published 



m 



Bolton et al.l (|2008l ). as well as the three disk-galaxy lenses found in it. 
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Fig. 2. — Full SDSS spectra of the selected gravitational lens candidates (from Princeton/MIT 
SDSS Spectroscopy). The background galaxy's emission lines are shown by marks in these spectra, 
and are presented in detail in Fig. [31 Information on the gravitational lens candidates is given in 
Table El 
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Fig. 3. — Emission lines of the background galaxies in the SDSS spectra for the selected gravitational 
lens candidates (from Princeton/MIT SDSS Spectroscopy). Although the SDSS spectra are in 
vacuum wavelengths, we indicate the emission lines in their usual air wavelengths for clarity. The 
full SDSS spectra are presented in Fig. [2l Information on the gravitational lens candidates is given 
in Table U 
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Fig. 3. — (continued) 
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and anomalous emission lines which are none of the above. Knowing the redshift of the foreground 
galaxy from the SDSS database, we can flag the emission lines belonging to the foreground galaxy as 
well as any skyline residuals. The remaining, anomalous emission lines are matched with each other 
to search for combinations of wavelength intervals which would correspond to the emission lines 
of a galaxy at a different redshift than the foregroung galaxy. For the selected spectra presenting 
evidence of anomalous emission lines from another galaxy along the line of sight, we inspect visually 
that the anomalous emission lines are not residuals from the subtraction of bright skylines, which 
were not flagged as such because they are too far from the central wavelength of the bright skyline. 
We eliminate in this way the false detections, and keep only the robust candidates. We select as 
lens candidates those spectra showing at least three emission lines among [O II] AA3727, H/3 4861, 
[O III] 4959, [O III] 5007 and Ha 6563, that belong to a higher redshift galaxy. 

In this pilot study, we looked only for bright lenses, that is, candidates we can easily follow- 
up with ground-based observations. Therefore, we decided to use a signal-to-noise (S/N) ratio 
per pixel selection to find emission lines, which is directly obtained from the ratio of the SDSS 
continuum-subtracted and noise spectra. This method is sufficient to detect bright emission lines, 
and saves the computation time that an integrated S/N ratio selection would require. The limits 
we used for the S/N ratio per pixel were, S/N > 5 for the first emission line peak and then > 3 
for the other emission line peaks (that were at the same redshift as the first emission line peak). 
This first step of the selection results in 20 spectra presenting a higher redshift galaxy aligned along 
the line-of-sight. 

We apply a second selection based on lensing probability: systems with a large strong lensing 
cross-section will have a higher probability for the background galaxy to be multiply imaged, and 
not only magnified. Studies of gravitational lens systems have found that isother mal mass models 



are a good approximation for the lensing mass distribution of early-type galaxies ([Koopmans et al 



20061 ). As we target massive disk-galaxy lenses, all our candidates are likely to have a non-negligible 
bulge component, which accounts for the first order lensing potential, and is well described by an 
isothermal sphere mass model. For a singular isothermal sphere (SIS) model, the strong lensing 



Table 2. Lens candidates 



Name 


R.A. (J2000) 


Decl. (J2000) 


Platc-MJD-Fibcr 


zp 


ZB 




J0812+5436 


08 12 13.30 


+54 36 50.9 


spSpec-53384-1871-114 


0.121 


0.358 


16.16 


J0903+5448 


09 03 15.62 


+54 48 56.4 


spSpec-51908-0450-226 


0.112 


0.689 


17.57 


J0942+6111 


09 42 49.08 


+61 11 15.5 


spSpec-51910-0486-155 


0.124 


0.716 


16.88 


J1150+1202 


11 50 19.56 


+12 02 57.3 


spSpec-53142-1609-141 


0.105 


0.532 


17.05 


J1200+4014 


12 00 46.90 


+40 14 00.0 


spSpcc-53401-1976-183 


0.116 


0.575 


16.78 


J1356+5615 


13 56 16.53 


+56 15 06.1 


spSpcc-52797-1323-531 


0.122 


0.600 


17.45 


J1455+5304 


14 55 43.55 


+53 04 41.2 


spSpec-52674-1164-270 


0.107 


0.635 


17.22 


J1625+2818 


16 25 51.95 


+28 18 21.4 


spSpec-52822-1408-417 


0.119 


0.349 


16.50 



^De Vaucouleurs model SDSS (AB) magnitude. 
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cross-section is asL = with ve = 4 it (a'^/c^) {D^s / Ds), where r^; is the Einstein radius, 

a is the velocity dispersion of the lens galaxy, Dls is the angular diameter distance between the 
lens and the source, which depends on the lens redshift, the source redshift and the cosmological 
parameters, and Dg is the source angular diameter distance. The lensing cross-section increases as 
a function of the lens galaxy mass and of the lens-source angular diameter distance, that is, of the 
redshift interval between the lens and the source. Unfortunately, we do not have measurements of 
o", which provides a prior on the mass for the candidate galaxies in our sample. However, we can 
establish a qualitative comparison of our lens candidates to the confirmed SLACS lenses, which were 
selected using a SIS model, to determine the highest lensing cross-section candidates. In addition 
to the redshift interval between the lens and the source, we can use as a comparison criterion the 
observed r-band magnitude from the SDSS database to obtain a rough estimate of the old stellar 
population luminosity and hence the mass of each galaxy (that is, the most massive galaxies are 
those with brighter r magnitudes). 

We select a final sample of eight gravitational lens candidates as those systems with the same 
range in r magnitude, foreground galaxy redshift zp and background galaxy redshift zb as the 
confirmed SLACS lenses (see Fig. [1]). The final sample has the following characteristics. All of the 
lens candidates have zp ~ 0.1. The background galaxy redshifts are distributed in two subsamples. 
The first in a group at 0.5 < zb < 0.8 with 16.7 < r < 17.7, and the second in a group of two lens 
candidates with specially low background galaxy redshifts at 0.3 < zb < 0.4, but more massive 
lens galaxies with 16 < r < 16.7. The list of final lens candidates with redshift and r magnitude 
information is presented in Table [2j The full SDSS spectra of the lens candidates are presented in 
Fig. [21 and the detected background emission lines are shown in Fig. [3l The twelve lens candidates 
which were discarded in the second selection step are shown in Appendix A. 



3. Follow-up of lens candidates 
3.1. Strategy 

We used a combination of u'-band imaging and spatially resolved long-slit spectroscopy, ob- 
served at the 2.5-m Nordic Optical Telescope (NOT), at the Observatorio del Roque de los Mucha- 
chos (Canary Islands, Spain), to confirm or reject the lensing hypothesis for the eight disk-lens 
candidates. A combination of both methods was required as the multiple images we want to detect 
are faint (they are not visible in the SDSS images) and are expected to be close to the galaxy 
center (disk galaxies have smaller Einstein radii than elliptical galaxies, due to their lower central 
concentration) . 

Imaging in the u' band provided, after modeling and subtraction of the foreground galaxy, an 
indication of the position of the potential lensed images, which guided us for positioning the slit for 
spectroscopic observations. However, the residuals in the image coul d also be non-smooth parts from 



the galaxy (an interesting discussion on that topic can be found in [Marshall et al.ll2008l ): therefore 
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we considered them only as potential lensed images of a background source, and upgraded them 
to confirmed Icnscd images only when finding in the long-slit spectroscopy a background galaxy 
emission line at the same position. The u' band was chosen over other optical bands because it 
minimizes the contamination from the foreground massive disk galaxy, which is faint in the blue 
end of the spectrum, while the background source is likely to be bright in the u' band, as indicated 
by its strong emission lines tracing the presence of star-forming regions. 

The SDSS spectra, which were taken with a 3" diameter fiber, have already shown that there 
is a background source near to the line-of-sight of the candidate lens galaxy. Spatially resolved 
long-slit spectroscopy allows us to detect spatially resolved multiple images of the emission lines of 
the background galaxy. As the slit can probe only a limited region around the foreground galaxy 
(e.g., the north-south axis, or the east-west axis), we placed it according to the galaxy-subtracted 
u'-band image. If we detected one potential lensed image we placed the slit in a position going 
through the image and the center of the galaxy; if we found multiple images around the galaxy, we 
took spectra in two different positions of the slit which optimized the number of potential lensed 
images covered by our observations. 

This strategy gave us a high chance of detecting the presence of multiple images of the same 
background source around the foreground lensing galaxy. 

Confirming the presence of a disk in the candidate lensing galaxies was achieved by taking 
additional optical images in the g', r' and i' bands. Four of the candidate lensing galaxies were 
visually identified as disk galaxies. For the other galaxies, we made a fit of the main galaxy 
component in the u' , g' , r' and i' bands and deduced a first-order estimate about the morphology 
of the galaxy. 

We now present the optical imaging and spectroscopy of the eight gravitational lens candidates. 



Table 3. Seeing (FWHM) conditions during observations (in arcsec) 



Name 


u' 


g' 


r' 


i' 


spectroscopy: 


spectroscopy: 




exptime 1800 s 


exptime 300 s 


exptime 300 s 


exptime 300 s 


east-west slit 


north-south slit 


J0812+5436 


1.15 


1.43 


1.39 


1.65 


1.08 




J0903+5448 


1.13 


1.3 


1.32 


1.50 




1.14 


J0942+6111 


1.30 


1.39 


1.24 


1.52 


0.70 




J1150+1202 


1.17 


2.02 


1.24 


1.71 


0.89 




J1200+4014 


1.17 


1.11 


1.17 


1.02 


0.76 


1.14 


J1356+5615 


1.11 


1.17 


1.15 


0.85 


0.76 


0.91 


J 1455+5304 


1.02 


1.28 


1.04 


0.95 




1.00 


J1625+2818 


1.00 








0.57 


0.86 
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3.2. Imaging 

3.2.1. Observations 

Optical images of the lens candidates were taken at the NOT with the MOSaic CAmera 
(MOSCA) during 2007 April 14-17 under variable weather conditions. We used the MOSCA in- 
strument to take advantage of its good sensitivity in the u' band. As the size of our targets does not 
require a mosaic camera, we centered our observations on one of the four CCDs of the mosaic, the 
CCDIO which showed the best properties regarding flux gradient and bad pixels on the chip. We 
used the SDSS filters u' , g' , r' and i' so as to be able to directly calibrate our observations with the 
SDSS database. The n'-band observations were performed during dark sky conditions and when the 
seeing conditions were favorable, in order to resolve the multiple images of the background sources. 
Exposure times of 1800 s were used due to the expected faintness of the lensed images. The g'-, r'- 
and i'-band observations were performed during grey sky conditions with exposure times of 300 s 
in each band. Seeing conditions during the observations are detailed in Table [3l the mean seeing 
during the imaging run was 1.26". 



3.2.2. Data reduction 



The optical images of the lens candidates were processed in IDL using standard data reduction 
techniques. The data were bia, s-subtracted and flat-fielded, the cosmic rays were removed using 
L. A. Cosmic (jvan Dokkumll200ll ). The dark current was negligible. 



3.2.3. Galaxy modeling 



The surfac e-brightness profi le of the galaxies were model ed using the publi cly available soft- 
ware GALFIT (jPeng et al.l 120021 '). We used a Sersic function (|Sersidll963l . Il968l ). that is, a seven 
parameter model (x, y, orientation, ellipticity, Sersic index n, scale radius and total magnitude), 
convolved with the point spread function, as defined from the stars in the field. The advantage of 
fitting a Sersic function is to recover a large range of galaxy profiles depending on the index n, for 
example, from the de Vaucouleurs profile (n = 4) to the exponential disk (n = 1). 

In the u' band, we use from one to three Sersic components in our models, depending on each 
case, in order to get a proper subtraction of the galaxy to identify the potential lensed images in the 
residuals. During this process, any strong emission from the candidate lensed background source 
was masked out to avoid skewing the fit to the lensing galaxy. The u'-band images and residual 
images are presented in Figs. [4] to [TT] (a and h). The results are summarized case by case in § 13.41 

The data for the lensing galaxies in the g' , r', and i'-bands were modeled by only one Sersic 
component, as here we are only interested in the morphology of the main component. This part is 
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Fig. 4. — Panel a) shows the w'-band image of J0812+5436, and panel b) the u'-band residual image after 
galaxy subtraction. The dashed lines indicate the position and width of the slit used to take long-slit spectra 

of the target. Panel c) shows a section of the 2D spectrum of .J0812+5436 with the [O II] AA3727 emission 
line of the background galaxy, at A = 5062 A, after skyline and galaxy subtraction, convolved with a gaussian 
of cr = 1 pixel for display purposes ; panel d) shows the same part of spectrum after skyline subtraction only. 
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Fig. 5. — Panel a) shows the w'-band image of J0903-I-5448, and panel b) the w'-band residual image after 
galaxy subtraction, convolved with a gaussian of rr = 2 pixels for display purposes. The dashed lines indicate 
the position and width of the slit used to take long-slit spectra of the target. Panel c) shows a section of the 
2D spectrum of J0903-I-5448 with the [O II] AA3727 emission line of the background galaxy, at A = 6295 A, 
after skyline and galaxy subtraction, convolved with a gaussian of cr = 2 pixels for display purposes; panel 
d) shows the same part of spectrum after skyline subtraction only. 
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Fig. 6. — Panel a) shows the u'-band image of J0942-I-6111, and panel b) the w'-band residual image after 
galaxy subtraction, convolved with a gaussian of ct = 2 pixels for display purposes. The dashed lines indicate 
the position and width of the slit used to take long-slit spectra of the target. Panel c) shows a section of the 
2D spectrum of J0942-I-6111 with the [O II] AA3727 emission line of the background galaxy, at A = 6397 A, 
after skyline and galaxy subtraction; panel d) shows the same part of spectrum after skyline subtraction 
only. 
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Fig. 7. — Panel a) shows the u'-band image of Jl 150+1202, and panel b) the w'-band residual image after 
galaxy subtraction, convolved with a gaussian of ct = 1 pixel for display purposes. The dashed lines indicate 
the position and width of the slit used to take long-slit spectra of the target. Panel c) shows a section of the 
2D spectrum of J1150+1202 with the [O 11] AA3727 emission line of the background galaxy, at A = 5710 A, 
after skyline and galaxy subtraction; panel d) shows the same part of spectrum after skyline subtraction 
only. 
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Fig. 8. — Panels a) and e) show the u'-band image of J1200-I-4014. Panels b) and f) show the w'-band 
residual image after galaxy subtraction, convolved with a gaussian of tr = 2 pixels for display purposes. 
The dashed lines indicate the position and width of the slit used to take respectively the east-west and 
north-south long-slit spectra of the target. Panels c) and g) show a section of, respectively, the east-west 
and north-south 2D spectra of J1200-I-4014, with the [O II] AA3727 emission line of the background galaxy 
at A = 5871 A. The spectra were skyline- and galaxy-subtracted, and convolved with a gaussian of a = 1 
pixel in panel c) and cr = 2 pixels in panel g) for display purposes. Panel d) and h) show the same part of 
the spectra after skyline subtraction only. 
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Fig. 9. — Panels a) and e) show the u'-band image of J1356-I-5615. Panels b) and f) show the w'-band 
residual image after galaxy subtraction, convolved with a gaussian of tr = 3 pixels for display purposes. 
The dashed lines indicate the position and width of the slit used to take respectively the east-west and 
north-south long-slit spectra of the target. Panels c) and g) show a section of, respectively, the east- west and 
north-south 2D spectra of J1356-I-5615, with the [O II] AA3727 emission line of the background galaxy at 
A = 5962 A. The spectra were skyline- and galaxy-subtracted, and the spectrum in panel c) was convolved 
with a gaussian of cr = 2 pixels for display purposes. Panel d) and h) show the same part of the spectra 
after skyline subtraction only. 
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Fig. 10. — Panel a) shows the u'-band image of J1455-I-5304, and panel b) the u'-band residual image after 

galaxy subtraction, convolved with a gaussian of ct = 1 pixel for display purposes. The dashed lines indicate 
the position and width of the slit used to take long-slit spectra of the target. Panel c) shows a section of the 
2D spectrum of J1455-I-5304 with the [O II] AA3727 emission line of the background galaxy, at A = 6095 A, 
after skyline and galaxy subtraction; panel d) shows the same part of spectrum after skyline subtraction 
only. 



- 21 - 




Fig. 11. — Panels a) and e) show the u'-band image of J1625+2818. Panels b) and f) show the u'-band 
residual image after galaxy subtraction, convolved with a gaussian of a = 2 pixels for display purposes. 
The dashed lines indicate the position and width of the slit used to take respectively the east-west and 
north-south long-slit spectra of the target. Panels c) and g) show a section of, respectively, the east- west 
and north-south 2D spectra of .J1625+2818, with the [O II] AA3727 emission line of the background galaxy 
at A = 5029 A. The spectra were skyline- and galaxy-subtracted, and convolved with a gaussian of a = 1 
pixel for display purposes. Panel d) and h) show the same part of the spectra after skyline subtraction only. 
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treated in more detail in § 13. 5[ 

3.3. Long-slit spectroscopy 

3.3.1. Observations 

Long-slit spectroscopy of the lens candidates was obtained at the NOT with the Andalucia 
Faint Object Spectrograph and Camera (ALFOSC) during 2007 May 7-11 under variable weather 
conditions. We used a slit width ranging from 1" to 1.3" depending on the seeing. We chose the 
grism #4 which covers a wavelength range 3200 A to 9100 A matching that of the SDSS spectra. 
This was a crucial point in order to confirm the presence of the three or more emission lines detected 
in the SDSS spectra of our gravitational lens candidates. However, it was at the expense of the 
resolving power of R~710, which is lower than that used for the SDSS spectra. The grism #4 has 
a blaze angle of — 1.5887 deg and a blaze wavelength of 496 nm. 

Seeing conditions during the observations are detailed in Table [3l the mean seeing during the 
spectroscopic run was 0.89". 

The slit position was determined based on the results of the galaxy subtraction presented in 
§ 13. 4[ For simplicity, we used only two slit positions (north-south and east-west) to cover the 
potential lensed images. 

3.3.2. Data reduction 

The spectroscopic data were processed using standard data reduction methods with the Image 
Reduction and Analysis Facility (IRAF), and using L. A. Cosmic for the cosmic-ray removal. The 2D 
spectra were bias subtracted, fiat-fielded, and dispersion corrected using He and Ne arc spectra. The 
skylines were subtracted by fitting a low order cubic spline interpolation to each row and subtracting 
it from the 2D spectrum. The galaxy contribution was subtracted using a similar technique along 
the columns of the galaxy spectrum. ID spectra were extracted from the calibrated, skyline- 
subtracted 2D spectra, in order to measure with better precision the wavelength of the background 
source emission lines, and the redshift of the background source. 

3.4. Combined results from imaging and long-slit spectroscopy 

We present in a combined figure, for each lens candidate, the u' band photometry and the 
long-slit spectroscopy, showing both the observations and the residuals after galaxy subtraction 
(see Figs. Hltollip. The optical images and the 2D spectra are presented with an identical spatial 
scale for the horizontal axis. The slit is indicated for each galaxy by two dashed lines in the u'- 
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band images. In the 2D spectra, the spatial axis corresponds to the position along the slit, and 
the wavelength axis corresponds to a wavelength range around the expected [O II] AA3727 emission 
line of the background galaxy. When spectroscopic observations of a lens candidate were obtained 
in two orientations of the slit, the two spectra are presented in the same figure, with the galaxy 
image rotated to present the slit horizontally, as indicated by the compass on each image. Below 
we discuss the lens candidates case by case, including the photometric data, spectroscopic data and 
combined results. 

3.4.1. J0812+5436 

Imaging-The galaxy-subtracted li'-band image shows two bright residuals linked by an arc 
around the north of the galaxy, suggestive of lensed images (see Fig. HI a and b). The two structures 
at the south of the galaxy are difficult to identify, but the nearest to the galaxy might be a lensed 
image. We planned spectroscopic observations in two positions of the slit: an orientation east-west 
passing through the side images, and an orientation north-south to study the southern objects. 
Unfortunately, bad weather conditions prevented us from taking the north-south spectrum. 

2D spectrum-We study the 2D spectrum around A = 5062 A, the wavelength of the expected 
[O II] AA3727 emission line. After galaxy subtraction, two bright residual images appear on each 
side of the galaxy at the wavelength A = 5060 A (see Fig.[31 c and d). The west image is at 0.7 ±1.1" 
and the east image at 1.3 it 0.6" from the center of the galaxy spectrum. They correspond to the 
positions through the slit of the images detected in the u' band. Although the two images in the 
spectrum are not neatly separated, the faint signal between them corresponds to the arc detected 
in the u' band, thus confirming the strong lens nature of the candidate. 

Conclusions-We confirmed that J0812+5436 is a strong gravitational lens. The Einstein radius 
is rE = 1.8 lb 0.4" (3.9 it 0.9 kpc), based on the half distance between the two confirmed lensed 
images, in the u' band image. The redshift of the lensed galaxy is zb = 0.357, as measured from the 
one dimensional spectrum. However, it is not possible to detail the nature of the residual images 
at the south of the galaxy, as we lack spectroscopic observations for them. The object nearest to 
the galaxy might be a lensed image of the same background source, while the furthest object might 
be unrelated. 

3.4.2. J0903+5448 

Imaging-The galaxy-subtracted u'-band image shows two residuals suggestive of lensed images 
(see Fig. \5\ a and b). The north image is at 0.6 it 0.6" and the south image at 0.9 it 0.3" from the 
galaxy center. We obtained spectroscopic observations with the slit in a north-south orientation, 
going through the two images. 
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2D spectrum-We study the 2D spectrum around A = 6295 A, the wavelength of the expected 
[O II] AA3727 emission Hue. After galaxy subtraction, a faint and continuous signal across the 
galaxy width appears at A = 6298 A (see Fig. [5l c and d), extending from 1.0" north to 2.0" south 
of the galaxy center. However, we do not detect distinctly an emission line from the background 
galaxy, nor if it is multiply imaged. The proximity of a strong skyline at A = 6300 A, the subtraction 
of which is bringing noise in this region of the spectrum, makes the study difficult. The [O III] 5007 
emission line would in theory bring better information on the lens candidate, presenting a better 
S/N ratio (see Fig. [3]), but in practice the fringing in the red part of our spectra makes it even 
more difficult to study. 

Conclusions-The presence of a background galaxy near J0903+5448 is confirmed from the 
SDSS (Princeton/MIT spectroscopy) one dimensional spectrum of the target, at a redshift zb = 
0.689. However, we do not detect properly the emission lines from the background galaxy in our 
own spectrum, due to the vicinity of a strong sky line near to the [O II] AA3727 emission line, and 
the fringing in the red part of the spectrum where the other strong emission lines of the background 
galaxy are found. Although the two residual images appearing in the u' band observation are very 
suggestive of lensing, better spectroscopic observations of the [O III] 5007 emission line would be 
necessary to confirm the lens nature of the candidate. 

3.4.3. JO942+6III 

Imaging-The galaxy-subtracted li'-band image shows a residual image at the east of the galaxy 
center (see Fig. [U a and b). We planned spectroscopic observations with the slit in the east- west 
orientation, passing through the potential lensed image and the center of the galaxy. 

2D spectrum-We study the 2D spectrum around A = 6397 A, the wavelength of the expected 
[O II] AA3727 emission line. After galaxy subtraction, one bright residual image appears at 1.1 ±0.8" 
east of the galaxy center, at the wavelength A = 6396 A (see Fig. [6l c and d). 

Conclusions-We confirmed the presence of an image of a background source near J0942+6111. 
The redshift of the background galaxy is zb = 0.716, as measured from the one dimensional 
spectrum. We can not confirm that the system is a lens, as we detect only one image of the 
background source. 

3.4.4. J1150+1202 

Imaging-The galaxy-subtracted u'-band image shows three bright residual images in a west, 
north-east and south-east configuration (see Fig. [71 a and b). We planned spectroscopic observations 
with east-west and north-south orientations of the slit to cover the three potential lensed images. 
Unfortunately, bad weather conditions prevented us from taking the north-south spectrum. 
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2D spectrum-We study the 2D spectrum around A = 5710 A, the wavelength of the expected 
[O II] AA3727 emission hne. After galaxy subtraction, one bright residual image appears at 0.8 ±0.8" 
west of the galaxy center, at the wavelength A = 5704 A (see Fig. [71 c and d). The emission line 
is superposed to the Mg I 5175 absorption line (in white in Fig. [TJ c) of the foreground galaxy at 
redshift Zf = 0.105. A faint extension of the [O II] AA3727 emission line is seen at the east of the 
galaxy, which may indicate a second image of the emission line. Comparing to the u' band, the 
bright image of the emission line corresponds to the west image, and the faint extension corresponds 
to the north-east image. As the west image appears brighter than the other image in the u' band 
observations, this explains that it is less affected by the presence of the Mg I 5175 absorption line. 

Conclusions-We confirmed the presence of an image of a background source near J1150+1202. 
The redshift of the background galaxy is zb = 0.531, as measured from the one dimensional 
spectrum. We can not confirm that the source is multiply imaged, but the presence of a faint 
extension to the [O II] AA3727 emission line concording with the north-east image in the n'-band 
residuals gives a strong indication for multiple images. It is not possible to detail the nature of the 
residual image at the south of the galaxy, as we lack spectroscopic observations for it. However, 
the presence of three images in the it'-band residual image, concording with one and possibly two 
spectroscopically confirmed images of the background source, while the third image is not proven 
spectroscopically, indicates the system is probably a strong gravitational lens. 



3.4.5. J 1200+40 14 

Imaging-The galaxy-subtracted n'-band image shows three residual images in a north-east, 
north-west and south configuration, suggestive of lensed images (see Fig. [8l a and b). Moreover, 
the three images seem to be linked by a faint ring, suggestive of an Einstein ring superposed 
to the lensed images. We obtained spectroscopic observations in the east-west and north-south 
orientations, covering the three potential lensed images. 

2D spectrum~We study the 2D spectra around A = 5871 A, the wavelength of the expected 
[O II] AA3727 emission line. After galaxy subtraction in the east-west spectrum, two bright residual 
images appear on each side of the galaxy at the wavelength A = 5870 A (see Fig. [HI c and d). In 
spite of the presence of a strong skyline at A = 5890 A, the subtraction of which induces noise in 
this part of the spectrum, the presence of the two residual images at the wavelength of the expected 
[O II] AA3727 emission line gives a strong indication for multiple images of the background source. 
The east image is at 1.7 ±0.7" and the west image at 1.0 ±0.6" from the center of the galaxy 
spectrum. Comparing the position of the two images to the ti'-band residual image, it appears 
these two images correspond to the north-east and north-west images. After galaxy subtraction in 
the north-south spectrum, we detect a faint image of the source on the north side of the galaxy at 
the wavelength A = 5870 A, but we find no evidence for the bright south image, which is therefore 
not belonging to the same source (see Fig. [Sj g and h). 
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Conclusions-We confirmed the presence of a background galaxy near J1200+4014, with a 
strong indication that it is multiply imaged, and therefore, that the system is a lens. The redshift 
of the background galaxy is zb = 0.575, as measured from the one dimensional spectrum. The 
third image, at the south of the galaxy, is not an image of the identified background source, but 
probably an unrelated object which is almost superposed to the foreground galaxy. 

3.4.6. J1356+5615 

Imaging-The galaxy-subtracted n'-band image shows two residual images at the east and south 
of the galaxy (see Fig. [9l a and b). There is also a central image which might be an artifact due to 
the difficulty of fitting the central pixels. We obtained spectroscopic observations in the east-west 
and north-south orientations of the slit, covering the three potential lensed images. 

2D spectrum-We study the 2D spectra around A = 5962 A, the wavelength of the expected 
[O II] AA3727 emission line. After galaxy subtraction of each of the east-west and north-south 
spectra, a residual image appears in each spectrum, respectively at 0.7 it 0.5" west and 0.7 ± 0.3" 
north of the center of the galaxy spectrum, at the wavelength A = 5962 A (see Fig. ^ c, d, g and 
h). These correspond to the central image in the w'-band residual image. We do not detect any of 
the two other images present in the u'-band residual image. 

Conclusions-We confirmed the presence of a central image of a background galaxy detected 
in both orientations of the slit. The redshift of the background galaxy is zb = 0.575, as measured 
from the one dimensional spectrum. It is not possible to determine if the images in the two spectra 
are two detections of the same image, or if there are actually two central images of the background 
source, or if they form a lensing arc. Therefore we can not conclude whether J1356+5615 is a 
strong lens or not. 

3.4.7. JI455+53O4 

Imaging-The galaxy-subtracted ti'-band image shows two residual images suggestive of lensed 
images, one on each side of the galaxy center (see Fig. 1101 ^ ^iiid b). The north image is at 0.5 ±0.5" 
and the south image at 1.0 it 0.4" from the galaxy center. We obtained spectroscopic observations 
with the slit in the north-south orientation, going through the two images. 

2D spectrum — We study the 2D spectrum around A = 6095 A, the wavelength of the expected 
[O II] AA3727 emission line. After galaxy subtraction, a residual image extending across the galaxy 
width appears at the wavelength A = 6090 A (see Fig. [lOl c and d). Although we observe a 
continuous signal, we can make the hypothesis that there are two unresolved images, and estimate 
their positions assuming they are separated at the center of the galaxy. This leads to a north 
image at 0.6 it 0.6" and a south image at 0.9 it 0.9" from the center of the galaxy spectrum. This 
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is consistent with the position through the sht of the two images detected in the u' band, thus 
confirming the strong lens nature of the candidate. 

Conclusion~We confirmed that J1455+5304 is a strong gravitational lens. The Einstein radius 
is rE = 0.7 lb 0.4" (1.4 it 0.8 kpc), based on the half distance between the two confirmed lensed 
images, in the u' band image. The redshift of the background galaxy is zb = 0.634, as measured 
from the one dimensional spectrum. 

3.4.8. .J 1625 +28 18 

Imaging-The galaxy-subtracted -u'-band image shows residuals suggestive of lensing: one bright 
image at the south-west of the bulge, and three smaller images superposed to what seems an arc 
on the north-east side of the bulge (see Fig. [TTl a and b). We obtained spectroscopic observations 
in the east-west and north-south orientations, covering the potential lensed images. 

2D spectrum-We study the 2D spectra around A = 5029 A, the wavelength of the expected 
[O II] AA3727 emission line. After galaxy subtraction in the east-west spectrum, a residual image 
appears at the wavelength A = 5032 A, at 0.8 it 0.3" west from the center of the galaxy spectrum 
(see Fig. [TTl c and d). It corresponds to the position of the south-west image in the u'-band 
residual image. After galaxy subtraction in the north-south spectrum, a residual image appears 
at the wavelength A = 5032 A, at 0.9 it 0.5" south from the center of the galaxy spectrum (see 
Fig. [TTj g and h). When compared to the u'-band residual image, the position of the image in the 
north-south spectrum does not correspond to the position of the north-east images which the slit 
was intended to probe, but corresponds to the position of the south-west image at the edge of the 
slit. However, we can see on the u' band image that the south-west image is brighter and larger 
than the other residual images. Moreover, the spectrum of the foreground galaxy shows many 
emission lines: this makes the peak of the [O II] AA3727 emission line of the background galaxy 
difficult to detect among the peaks of the foreground galaxy (see Fig. [3]). This could explain that 
we detect only the brightest of the images of the background emission line, while the fainter images 
would be at the level of the noise occasioned by the foreground galaxy. 

Conclusions-We confirmed the presence of an image of a background source near J1625+2818. 
The redshift of the background galaxy is zb = 0.350, as measured from the one dimensional 
spectrum. We can not confirm that the source is multiply imaged, as we detect spectroscopically 
only the brightest image. The remaining residual images might be too faint to be detected in our 
spectra. However, the configuration of the residuals in the u' band image points toward a good 
probability that the system is a strong gravitational lens. 




Fig. 12. — Imaging in the r' band of the targets which arc obvious disk galaxies (J1625+2818 is 
shown in the u' band). The size of the images is 21.7" x 21.7". 
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Fig. 13. — Imaging in the r' band of the targets whose morphology can not be deduced by visual 
inspection and require galaxy modeling to know if they are disk galaxies. The indications of the 
galaxy morphology refer to the conclusions of § 13. 51 The size of the images is 21.7" x 21.7". 
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3.5. Morphology 

The galaxies J0903+5448, J0942+6111 and J1356+5615 can be identified as disk galaxies 
directly from looking at their optical imagej^. As for the galaxy J1625+2818, we can see it is a 
barred spiral with a dominant bulge (see Fig. [T2|) . However, the candidate massive disk galaxies 
J0812+5436, J1150+1202, J1200+4014 and J1455+5304 do not have characteristics that allow us 
to identify their galaxy type by simple visual inspection (see Fig. [T3]l . 

To gain further insight into the morphology of these galaxies, we used GALFIT to model the 
main component of each galaxy with a Sersic profile. We fitted the galaxies in each of the available 
photometric bands: u' , 5', r' and i' . For each fit we used a simple method which allows us to get 
robust information about the morphology of the main component of the galaxies in each band: we 
fitted the data with a Sersic profile keeping n fixed to 4 (de Vaucouleurs bulge), then fitted the 
data with a Sersic profile keeping n fixed to 1 (exponential disk). The comparison of the of the 
two fits tells us if the morphology of the galaxy is nearer to a bulge or a disk. We also fitted the 
data with a Sersic profile where n is let free. 

While this method allows us to determine if the galaxy morphology is definitely a disk or a 
bulge, it is difficult to identify SO galaxies, which have a Sersic index between the two models. For 
most of the cases, we can only say that the galaxy is not dominated by a disk, but we can not 
determine clearly between SO and elliptical galaxies. It is to note that even with HST imaging, it 
can be difficult to differentiate SO galaxies from elliptical galaxies, if the disk is not a dominant 
feature. 

Here follows a summary galaxy by galaxy of the results, which are presented in Table HI 
J 08 12+ 5436 

In the u' and g' bands the of the fit is smaller for a disk (n = 1) while in the r' and i' band the 
of the fit is smaller for a bulge (n = 4) . This indicates that the galaxy is composed of a disk and 
a bulge. However, due to the very bright lensed images superposed to the galaxy, the fit may not 
be completely reliable, although the lensed images were masked. Therefore, we can only conclude 
that J0812+5436 is a probable disk or SO galaxy. 

J 1150+ 1202 

In the u' band the of the fit is similar for a disk (n = 1) and a bulge (n = 4), while in the g\ r' 
and i' bands the of the fit is smaller for a bulge (n = 4). In addition, the best fit values in all 
the bands are for 1 < n < 4, which is between a disk and a bulge model. We conclude that SDSS 
J1150+1202 is probably an SO galaxy. 

J 1200+40 14 



Fitting these galaxies would not bring more information about their morphology, as GALFIT does not include 
the inclination of the disk as a fitting parameter. Therefore, the Sersic index n would be increased artificially, as the 
light received is more concentrated due to the inclination of the disk. 



- 31 - 



In all bands the of the fit is smaller for a bulge (n = 4). We conclude that SDSS J1200+4014 
is probably an SO or an elliptical galaxy. 

J 1455+ 5304 

In all bands the of the fit is smaller for a bulge (n = 4). We conclude that SDSS J1455+5304 
is probably an SO or an elliptical galaxy. 



4. Mass-to-light ratios 

We compute estimates of the total M/L ratios for the gravitational lenses we have discovered. 
For any circular lens, the lens potential is a function only of the distance from the lens center. The 
mass enclosed within the Einstein radius r^; is obtained from the relation r% 



The luminosity inside the Einstein radius is computed using a Sersic profile (iGraham &: Driver 



2OO5I ). O ur images are cal i brated using the SDSS database magnitudes, and K-correction is applied 



following lOke &: Sandagd (|l968l ). The results are presented in Table [5l However, these values are 
difficult to compare to usual M/L ratios, as they probe the total mass and not just the stellar or 
baryonic matter. 

In the case of the SO or elliptical galaxy lens J1455+5304, the Einstein radius is small {rE = 
0.7 ±0.4", that is, 1.4 ±0.8 kpc), so we can expect that the mass probed within is largely dominated 
by b aryons. On this basis, t he Mtot/Lj = 1.5 ± 0.9 for J1455±5304 is consistent with the value 



from lTrott fc Websteil pOOd ) of M/Lj = 1.1 ± 0.2 for the disk of the Sab spiral lens the Einstein 
Cross. 

In the case of the (probable disk-galaxy) lens J0812±5436, we obtain a particularly high M/L 
ratio, with M^ot/Li = 5.4 ± 1.5. The Einstein radius is larger than expected for galaxy-galaxy 
lensing systems, with ve = 1-8 ± 0.4" (3.9 ± 0.9 kpc). The total mass within might contain a 
large amount of DM, although it is still in the few central kiloparsecs of the galaxy. We note that 
J0812±5436 is surrounded by five galaxies at a similar redshift, with a group radius of ~ 200 kpc 
(see Fig. [TH and Table [6]) . This group of galaxies may add to the gravitational potential deflecting 
the background source. 

The limited observational precision that we obtain on measuring the Einstein radius and the 
surface brightness of the different components of the galaxies limits us in the study of the M/L 
ratios of these new galaxy-galaxy lenses. Therefore, we can not gain insight on the M/L ratios of 
the different components, nor constrain the shape of the IMF or the maximum disk hypothesis. 
HST or Adaptive Optics imaging is required to further study the M/L ratio of these lens galaxies. 
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Table 4. Morphology 







u' 


o' 


r' 


i' 








J0812+5436 






xl for n = 


1 


1.548 


1.378 


1.689 


2.424 


xl for n = 


4 


2.227 


1.417 


1.673 


1.482 


best fit n 




0.74 


1.82 


1.92 


2.35 








J1150+1202 






xl for n = 


1 


1.129 


1.129 


1.499 


1.279 


xl for n = 


4 


1.127 


1.110 


1.367 


1.172 


best fit n 




1.98 


2.32 


2.42 


3.06 








J1200+4014 






xl for n = 


1 


1.161 


1.453 


1.664 


1.697 


xl for n = 


4 


1.133 


1.204 


1.265 


1.264 


best fit n 




3.38 


2.83 


2.68 


2.82 








J1455+5304 






for n = 


1 


1.181 


1.216 


1.419 


1.517 


xl for n = 


4 


1.160 


1.162 


1.235 


1.357 


best fit n 




2.62 


2.35 


2.80 


2.52 



Table 5. Mass-to-light ratio within ve 



Name 


Mtot/Li 






i/Lb 


Mtot/Lu 


r-E 


(kpc) 


J0812+5436 


5.4 


± 1.5 


8.2 


±2.0 


11.2 


± 2.8 


14.6 


± 2.3 


3.9 


±0.9 


J1455+5304 


1.5 


±0.9 


2.2 


±1.3 


3.0 


±1.7 


3.1 


± 1.9 


1.4 


±0.8 



Table 6. Local galaxies around J0812+5436 (z=0.12P ) 



Name 


Photo z (SDSS) 


Projected distance from J0812+5436 


Galaxy A 


0.166 ±0.042 


~ 184 kpc 


Galaxy B 


0.135 ± 0.030 


~ 129 kpc 


Galaxy C 


0.133 ±0.010 


~ 158 kpc 


Galaxy D 


0.096 ± 0.008 


~ 219 kpc 


Galaxy E 


0.131 ±0.006 


~ 297 kpc 



''Spectroscopic redshift from SDSS. 
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5. Discussion 

The results of the fohow-up of the eight disk-galaxy gravitational lens candidates are summa- 
rized in Table [71 From the optical imaging and spectroscopy presented here, we can determine 
clearly that two systems are strong gravitational lenses, J0812-I-5436 being a probable disk or SO 
galaxy, and J1455-I-5304 being a SO or elliptical galaxy. 

Concerning the other candidates, we classify as probable gravitational lenses those systems 
in which we find indication of multiple imaging from both imaging and spectroscopy, although 
observations do not allow us a complete confirmation of the lens nature of the system. We find two 
probable gravitational lenses, J1150-I-1202 being probably a SO galaxy, and J1200-I-4014 being a SO 
or elliptical galaxy. 

Finally, we classify as unknown the systems for which we found indication of multiple imaging 
in optical observations, but could not confirm the presence of multiple images in the spectroscopy. 
These systems are good lens candidates but the quality of our observations did not allow us to 
conclude on their nature. In this category, we find J0903-I-5448, J0942-I-6111, and J1356-I-5615 
which are disk dominated galaxies, as well as the spiral galaxy J1625-I-2818. 

We see that for most of the gravitational lens candidates, our combined photometric and 
spectroscopic observations are not sufficient to determine clearly the presence of multiple images 
of the background source. This is due mainly to the fact that the images of the background source 
are very near to the center of the foreground galaxy, the Einstein radius being expected to be of 
the order of ~ 1". Therefore, ground-based observations suffer from seeing conditions for such a 
study. Moreover, some of the images might be obscured by dust in disk galaxies. 

Confirmation of the remaining lens candidates will require either high resolution spatially- 
resolved spectroscopy, or high resolution imaging. Integral field spectroscopy would allow us to 
probe both the presence of multiple images of a background source and the configuration of the 
lensed images. However, because the galaxies in our study are candidate massive disk galaxies, 
they are not expected to present much structures in the blue part of the spectrum, which would be 
mistaken for lensed images. In the case of the spiral lens J1625-I-2818, the position of the potential 



Table 7. Results 



Name 


Type 


Morphology 


J0812+5436 


strong lens 


probably disk 


J0903+5448 


unknown 


disk 


J0942+6111 


unknown 


disk 


J1150+1202 


probably strong lens 


probably SO 


J1200+4014 


probably strong lens 


SO or elliptical 


J1356+5615 


unknown 


disk 


J1455+5304 


strong lens 


SO or elliptical 


J1625+2818 


unknown 


spiral 
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Fig. 14. — Imaging in the r' band of the field around J0812+5436. The five circled galaxies have 
SDSS photometric redshifts close to the redshift of J0812+5436 (see Table [6|). The presence of a 
galaxy group may account for the large Einstein radius of J0812+5436 {te = 1-8 ± 0.4"). At the 
redshift z = 0.121, 30" correspond to 65 kpc. 



lensed images around the bulge also prevente d any contaminatioii by the arms for detecting them 



Ther efore HST imaging, or Adaptive Optics (jMarshall et al.ll2007l : iMcKean et al.ll2007l : ISluse et al 



20081 ) imaging, might be the best choice for confirming the genuine lens nature of our disk-galaxy 
lens candidates, confirming the disk nature of our new gravitational lenses, and studying the mass 
distribution of these systems. 

In the following we present suggestions for improving the efficiency of a search and for enlarging 
the size of the survey. 

1- Due to the very large number of galaxies in the SDSS Main Galaxy Sample, we covered only 
a part of the color parameter space corresponding to disk galaxies. We chose to limit ourselves to 
the region containing early- type disk galaxies, which are massive disk galaxies and therefore better 
suit our purpose. A larger survey should span all of the disk galaxy color range (with applying a 
suitable r magnitude criterion to select massive galaxies), although this would largely increase the 
size of the initial galaxy sample. 
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2- The automated method we used to select spectra showing evidence for a background galaxy 
along the line of sight is based on a S/N ratio per pixel emission-line detection. This was done in 
order to select bright peaks in the spectra. A complete search should be based on an integrated 
S/N ratio selection to detect and identify all background emission lines. 

3- We did not cover the redshift range z < 0.1, in order to preserve the efficiency of the search. 
Indeed, while half of the galaxies in the SDSS Main Sample lie in this region, very few lenses are 
expected at such a low redshift. However, it would be of interest for a larger survey to explore this 
redshift range. 

4- Selecting lens candidates based on the presence of three emission lines limits, in the case of the 
SDSS, the redshift of the background galaxy to z < 0.8. Selection based on one detected einission 
line and its identification as the [O II] AA3727 doublet, as in the OLS-survey (jWillis et al.ll2005l . 



20061 ). would allow us to enlarge the sample of lens candidates. Indeed, almost all our gravitational 



lens candidates have a resolved [O II] AA3727 doublet. However, it would also increase the false 
detection rate, hence reduce the efficiency of the survey. 

5- Among the SLACS lenses, some late- type galaxy lenses were found (IBolton et al.l 120081 ). al- 
though from targeting luminous red galaxies with absorption-line dominated spectra, expected to 
be mainly early- type galaxies. This can happen if the bulge in the late- type galaxy is prominent, 
and therefore this should be taken into account when searching for the disk galaxy population in 
other large surveys. 



6. Conclusion 



We have presented the first automated search for disk-galaxy gravitational lenses using the 
SDSS database, and proved the feasibility of such a project with the discovery of a galaxy-galaxy lens 
which is very probably a disk, of four interesting disk-galaxy lens candidates, and three confirmed 
or probable lenses which may be SO galaxies. This project is the first step in finding more disk- 
galaxy lenses in the SDSS, and possibly in other large surveys. To date, only ten confirmed 
disk-galaxy lenses are known. Assembling a larger sample of disk-galaxy lenses would open up 
promising perspectives for measuring the M/L ratios of disks and bulges, and for studying the 
structure of disk galaxies. HST or Adaptive Optics imaging will be needed to further study the 
systems presented in this work, particularly to probe the presence of a disk in the lens galaxy of 
J0812+5436 and study its mass distribution, as well as to confirm the lens nature of the interesting 
disk-galaxy lens candidates we discovered. 
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A. Additional disk-galaxy lens candidates 

During the spectroscopic selection of the disk-galaxy lens candidates (§ [2]), we found twenty 
systems with evidence of a background galaxy along the line-of-sight. Among these twenty disk- 
galaxy lens candidates, we found eight systems to have a similar distribution of redshifts and r 
magnitude as the SLAGS confirmed lenses (see Fig. [1]), and selected these eight systems for follow- 
up. We present in Appendix A the remaining twelve candidates: their redshifts and r magnitude 
distribution compared to that of the SLAGS lenses is shown in Fig. 1151 and details on the candidates 
are presented in Table El These candidates were discarded because they had both a short redshift 
interval between the background and the foreground galaxies, and a relatively faint r magnitude 
for the foreground galaxy. In comparison, the SLAGS lenses with such a short redshift interval 
show brighter r magnitudes (with the r magnitude being a rough tracer of the mass of the galaxy), 
which increases their lensing cross-section. 
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Fig. 15. — Distribution of the foreground galaxy redshift zp, the background galaxy redshift zb 
and the r magnitude of the twelve lens candidates we did not selected for foll ow-up, compared to 
the SLACS confirmed strong lenses. The sample of SLACS strong lenses from iBolton et al.l (j200(j| ) 
was the one we actually used to select our fina l lens candidates. We show also the final sample of 
SLACS strong lenses published in lBolton et al.l (|2008l ). as well as the three disk-galaxy lenses found 
in it. 



Table 8. Additional lens candidates 



Name Platc-MJD-Fiber zp zb r'^ 



SDSS 


J005621 


66-091201.9 


spSpec- 


52146- 


0658- 


393 





103 





313 


17 


21 


SDSS 


J090146 


27+554102.0 


spSpec- 


51908- 


0450- 


388 





116 





247 


17 


30 


SDSS 


J103143 


93+421859.3 


spSpec- 


53033- 


1360- 


415 





119 





308 


17 


49 


SDSS 


J110039 


10+120159.1 


spSpec- 


53119- 


1603- 


440 





129 





220 


17 


23 


SDSS 


J112500 


00+053604.7 


spSpec- 


52376- 


0836- 


536 





155 





319 


17 


60 


SDSS 


J114440 


12+043650.5 


spSpec- 


52373- 


0839- 


230 





104 





256 


17 


51 


SDSS 


J134308 


25+602754.8 


spSpec- 


52319- 


0786- 


236 





120 





320 


16 


86 


SDSS 


J145555 


90+040745.3 


spSpec- 


52045- 


0588- 


553 





126 





224 


17 


39 


SDSS 


J150339 


92+370728.1 


spSpec- 


52819- 


1352- 


472 





116 





369 


16 


89 


SDSS 


J171605 


16+275206.3 


spSpec- 


52410- 


0977- 


093 





119 





309 


17 


18 


SDSS 


J172614 


51+604142.5 


spSpec- 


51792- 


0354- 


151 





113 





230 


17 


28 


SDSS 


J235726 


45-090917.7 


spSpoc- 


52143- 


0650- 


370 





120 





394 


17 


45 



''De Vaucouleurs model SDSS (AB) magnitude. 



- 38 - 



REFERENCES 

Adelman-McCarthy, J. K., et al. 2007, ApJS, 172, 634 
Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212 
Belokurov, V., et al. 2007, ApJ, 671, L9 

Bolton, A. S., Buries, S., Schlegel, D. J., Eisenstein, D. J., & Brinkmann, J. 2004, AJ, 127, 1860 
Bolton, A. S., Buries, S., Koopmans, L. V. E., Treu, T., & Moustakas, L. A. 2006, ApJ, 638, 703 
Bolton, A. S., & Buries, S. 2007, New Journal of Physics, 9, 443 

Bolton, A. S., Buries, S., Koopmans, L. V. E., Treu, T., Gavazzi, R., Moustakas, L. A., Wayth, R., 
& Schlegel, D. J. 2008, ApJ, 682, 964 

Browne, I. W. A., et al. 2003, MNRAS, 341, 13 

Castander, F. J., Treister, E., Maza, J., & Gawiser, E. 2006, ApJ, 652, 955 
Covone, G., et al. 2008. rarXiv:0809.4125l 

de Jong, R. S., & Bell, E. F. 2007, Island Universes - Structure and Evolution of Disk Galaxies, 
107 

Dobler, G., Keeton, C. R., Bolton, A. S., & Buries, S. 2008, ApJ, 685, 57 
Dutton, A. A., van den Bosch, F. C., Dekel, A., & Courteau, S. 2007, ApJ, 654, 27 
Fassnacht, C. D., et al. 2006, ApJ, 651, 667 
Faure, C., et al. 2008, ApJS, 176, 19 

Flynn, C., Holmberg, J., Portinari, L., Fuchs, B., & Jahreifi, H. 2006, MNRAS, 372, 1149 
Fukugita, M., Shimasaku, K., & Ichikawa, T. 1995, PASP, 107, 945 
Gentile, G., Salucci, P., Klein, U., Vergani, D., & Kalberla, P. 2004, MNRAS, 351, 903 
Ghosh, K. K, & Narasimha, D. 2008, larXiv:0810.3971l 

Governato, F., Willman, B., Mayer, L., Brooks, A., Stinson, G., Valenzuela, O., Wadsley, J., & 
Quinn, T. 2007, MNRAS, 374, 1479 

Graham, A. W., & Driver, S. P. 2005, Publications of the Astronomical Society of Australia, 22, 
118 

Huchra, J., Gorenstein, M., Kent, S., Shapiro, I., Smith, G., Horine, E., &: Perley, R. 1985, AJ, 90, 
691 



- 39 - 



Inada, N., et al. 2008. larXiv:0809.0912l 

Jaunsen, A. O., & Hjorth, J. 1997, A&A, 317, L39 

Kassin, S. A., de Jong, R. S., & Weiner, B. J. 2006, ApJ, 643, 804 

Kayo, I., et al. 2007, AJ, 134, 1515 

Keeton, C. R., & Kochanek, C. S. 1998, ApJ, 495, 157 

Koopmans, L. V. E., de Bruyn, A. G., & Jackson, N. 1998, MNRAS, 295, 534 

Koopmans L. V. E., Treu T., Bolton A. S., Buries S., Moustakas L. A., 2006, ApJ, 649, 599 

Limousin, M., et al. 2008, larXiv:0812,1033l 

Mailer, A. H., Simard, L., Guhathakurta, P., Hjorth, J., Jaunsen, A. O., Flores, R. A., & Primack, 
J. R. 2000, ApJ, 533, 194 

Mandelbaum, R., van de Ven, G., & Keeton, C. R. 2008, arXiv:0808.2497] 

Marshall, P. J., et al. 2007, ApJ, 671, 1196 

MarshaU, P. J., Hogg, D. W., Moustakas, L. A., Fassnacht, C. D., Bradac, M., Schrabback, T., k 
Blandford, R. D. 2008. larXiv:0805.14 69 l 

McKean, J. P., et al. 2007, MNRAS, 378, 109 

Moller, O., Kitzbichler, M., & Natarajan, P. 2007, MNRAS, 379, 1195 

Myers, S. T., et al. 2003, MNRAS, 341, 1 

Navarro, J. F., & Steinmetz, M. 2000a, ApJ, 538, 477 

Navarro, J. F., & Steinmetz, M. 2000b, ApJ, 528, 607 

Newton, E. R., Marshall, P. J., & Treu, T. 2008. larXiv:0810.3934l 

Oguri, M., et al. 2006, AJ, 132, 999 

Oguri, M., et al. 2008, AJ, 135, 520 

Oguri, M. 2007, New Journal of Physics, 9, 442 

Oguri, M., Inada, N., Blackburne, J. A., Shin, M.-S., Kayo, I., Strauss, M. A., Schneider, D. P., & 
York, D. G. 2008, MNRAS, 1321 

Oke, J. B., & Sandage, A. 1968, ApJ, 154, 21 



-40- 



Patnaik, A. R., Browne, I. W. A., King, L. J., Muxlow, T. W. B., Walsh, D., & Wilkinson, P. N. 
1992, Gravitational Lenses, 406, 140 

Peng, C. Y., Ho, L. C, Impey, C. D., k Rix, H.-W. 2002, AJ, 124, 266 

Percival, W. J., et al. 2007, ApJ, 657, 645 

Portinari, L., Sommer-Larsen, J., & Tantalo, R. 2004, MNRAS, 347, 691 

Portinari, L., & Sommer-Larsen, J. 2007, MNRAS, 375, 913 

Salucci, P., & Persic, M. 1999, A&A, 351, 442 

Salucci, P., Yegorova, 1. A., & Drory, N. 2008, MNRAS, 388, 159 

Sersic, J. L. 1963, Boletin de la Asociacion Argentina de Astronomia La Plata Argentina, 6, 41 

Sersic, J. L. 1968, Cordoba, Argentina: Observatorio Astronomico, 1968 

Shin, E. M., & Evans, N. W. 2007, MNRAS, 374, 1427 

Sluse, D., Courbin, F., Eigenbrod, A., & Meylan, G. 2008, A&A, 492, 39 

Strateva, I., et al. 2001, AJ, 122, 1861 

Trott, C. M., & Webster, R. L. 2002, MNRAS, 334, 621 

Trott, C. M., Treu, T., Koopmans, L. V. E., & Webster, R. L. 2008. larXiv:0812.0748l 

Tully, R. B., & Fisher, J. R. 1977, A&A, 54, 661 

Vallejo, O., Braine, J., & Baudry, A. 2002, A&A, 387, 429 

van Albada, T. S., & Sancisi, R. 1986, Royal Society of London Philosophical Transactions Series 
A, 320, 447 

van den Bosch, F. C., Mo, H. J., & Yang, X. 2003, MNRAS, 345, 923 
van Dokkum, P. G. 2001, PASP, 113, 1420 

Willis, J. P., Hewett, P. C., & Warren, S. J. 2005, MNRAS, 363, 1369 

Willis, J. P., Hewett, P. C., Warren, S. J., Dye, S., & Maddox, N. 2006, MNRAS, 369, 1521 

Winn, J. N., Kochanek, C. S., McLeod, B. A., Falco, E. E., Impey, C. D., & Rix, H.-W. 2002, ApJ, 
575, 103 

Winn, J. N., Hall, P. B., k Schechter, P. L. 2003, ApJ, 597, 672 
York, D. G., et al. 2000, AJ, 120, 1579 



-41 - 

York, T., Jackson, N., Browne, I. W. A., Wucknitz, O., & Skelton, J. E. 2005, MNRAS, 357, 124 



This preprint was prepared with the AAS IM^^X macros v5.2. 



